function E=model_err(B)

%Recibe una matriz B que representa los pixeles (en tonos de gris) de una
%imagen y devulve la diferencia entre la prediccion del modelo y la imagen
%real E.

[M,N]=size(B);
aux_im=zeros(M+1,N+2);
aux_im(2:end,2:end-1)=B;
E=zeros(size(B));
for i=1:M
    line_ant=aux_im(i,:);
   for j=1:N
       a=aux_im(i+1,j);
       b=line_ant(j+1);
       c=line_ant(j);
       d=line_ant(j+2);
       E(i,j)=double(B(i,j))-double(predictorMED(a,b,c,d));
       if (a==b)&&(b==c)&&(c==d)&&(E(i,j)==0)
           E(i,j)=256;
       end
   end
end